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Abstract 

Abel Inversion is currently used in laser-plasma studies in order to estimate the 
electronic density Ue from the phase-shift map Sep obtained via interferometry. The 
main limitation of the Abel method is due to the assumption of axial symmetry of 
the electronic density, which is often hardly fulfilled. In this paper we present an 
improvement to the Abel inversion technique in which the axial symmetry condition is 
relaxed by means of a truncated Legendre Polinomial expansion in the azimutal angle. 
With the help of simulated interferograms, we will show that the generalized Abel 
inversion generates accurate densities maps when applied to non axisymmetric density 
sources. 

Abel Inversion is widely used in many context and, in laser-plasma studies it leads to 
a 2D electronic density map Ug reconstruction from phase-shift maps 5(j) recorded using 
interferometry Once the phase map 6(j){z, x) has been extracted from the interferogram 

via standard FFT techniques P] or with a more sophisticated Wavelet-Based method [Q, 
the best symmetry axis zq should be defined and two half-space phase maps 

6(f)'^{(,x) = 6(f>{z — Zo,x) Z > Zq 

6(f)~{(,x) = 6(j){Zo — Z,x) Z < Zq (0.1) 
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are introduced. By assuming axial symmetry around the laser-beam propagation axis x, 
50"*" = 50~ so the electronic density map ne(r, x) is computed as 



ne{r,x) 



= -n, 



1 



d 



(0.2) 



where Xp is the probe wavelength and ric = nm{c/{eXp))'^ is the critical density at A- 



Though laser plasmas show approximately axial symmetry in general, significant devi- 
ations from the symmetry may occur. In these latter cases the Abel Inversion applied to 
an artificial profile obtained by symmetrization of either 5(j)'^ or 50", can lead to mislead- 
ing reconstruction of the density distribution. Alternately, to consider the "mean" phase 
distribution = 1/2(50+ + 50") can also induce large errors. 

In 1981 Yasutomo et al ^ presented a generalization of Abel inversion based on the 
assumption that the 3D density distribution n{r,x, 6) can be factorized as the product of a 
isotropic density nQ{r) and a corrective term g{rsmd), being 9 the azimutal angle. 

In this paper we introduce a generalization of the Abel inversion algorithm to be applied 
to moderately asymmetric interferograms. Unlike Yasutomo, we base our algorithm on a 
Legendre polinomial expansion of n(r, x, 9) in the angular variable alone. We will show 
that such an extension of the Abel method allows accurate reconstructions of the density 
distribution in some simulated non symmetric cases. 

The basic geometry of the phase-shift acquisition via interferometry is shown in Fig. |1|, 
in which a plane parallel to the laser propagation axis (x-direction). Fig. |l] a), and a plane 
perpendicular to that axis. Fig. |I|b), are shown, respectively. 

Let us point out that any departure from the mirror symmetry respect to the plane 
perpendicular to the probe axis (the x — z plane in Fig. |ip cannot be taken into account 
because of the line integral in the acquisition step. Consequently we can assume that such 
a mirror symmetry is satisfied. 

The phase-shift 5(j){z, x) detected in the (z, x) position on the interferogram is then linked 
to the electronic density ne{r,9,x) as: 



where we moved to the cylindrical coordinates (r, 9), taking x as symmetry axis. 

We now make the physically justified assumption that the angular dependence of ne(r, 9, x) 
is everywhere "well behaved" (no abrupt changes occur) so that Ue can be developed as a 




(0.3) 
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truncated series of orthonormal Legendre Polinomials Pi{cos{6)): 



L 

ne(r, 9,x)=Y, Mr, x)Pi{cos{9)) . (0.4) 

1=0 

To find the appropriate value of L, let us simply observe that for each x from S(j)^{() and 
6(j)^{Q we can build up two linearly independent sequences 

50,^i(50+(C) + 5r(C)) 

5(/..^^(50+(C)-50-(C) 

so that for each x and r we have two independent degrees of freedom which could be linked 
to the angular dependence and this leads to L = 1. The truncation of the series of Eq. p.4| 
up to Z = 1 is straightforward. Since Po{x) = 1 and Pi{x) = x, we have 

n(r, 9, x) = no{r, x) + ni(r, x) cos{9) (0.5) 

and so the phase-shift map is computed as: 

, 2n r if,, , \{z ~ zq)\ ^ , . 

z,x) = —— I I — \nQ[r,x) + ni[r,x) ) dr . (0.6) 



P J\z-zo\ — (2; — z^y 

Defining Vi (r, x) = ni (r, x)/r and extracting the symmetric and antisymmetric components 
of 6(j) in Eq. p.6| , we obtain 

S(t)s{C,x) = -— / , ., ... dr 



7r(^ r v\{r^x) 



and we can finally invert Eq. p.7| obtaining the coefficients of the generalized Abel Inversion: 

X f°° Id 



The application of the generalized Abel inversion (equations ( |0]^), (p.8|)) is straightfor- 
ward and very effective. In order to prove this, the new algorithm will be tested with two 
sample interferograms. Both of them have been obtained numerically from a priori known 
density distributions, with which the reconstructed distributions can be compared. 
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Let us firstly consider the sample interferogram of Fig. ^ b), obtained from the 3D 
density distribution of the form n"^™'^(r, 6*, x) = no(r, x) Po(cos(^)) + ni{r,x) Pi{cos{9)) + 
n2{r,x) P2{cos{9)), with rii gaussian shaped in the radial direction r and exponentially de- 
creasing in the longitudinal direction x. The simulation is performed assuming a l/zm 
wavelength probe and a maximum density of nj in Uc units as n"^"-^ = 0.1, n^"^ = 0.05, 

^max ^ Q Q25_ 

Once the phase shift S(j) has been extracted from the interferogram, an automatic pro- 
cedure to optimize the position of the global symmetry axis (2:0) has been applied and the 
two half phase maps 50"'" and S(j)~ have been constructed. Next, the electronic densities 
n~^{r,x), n~{r,x), n"*^'*"(r, x) are computed via standard Abel Inversion applied to 6(f)~^, 
6(j)~ , 6(f)s = 1/2(50+ + 6(f)~), respectively. Finally, the generalized Abel inversion (Eqq. p75| 
p.8| ) is applied to both 6(j)^ and 6(f)~ producing n^'^^'^^{r,6,x). 

In Fig. H projections onto the z — x plane of the simulated density (a), the standard 
Abel inversion of the symmetryzed map Scps (c) and of the generalized Abel inversion of 5(j)'^ 
and 6(j)~ (d), are shown. As it is clear, n"^'^"'^ shape differs considerably from the one of 
j^True^ while n'^^"'^^ well match the true density map. For a more quantitative comparison, 
we have reproduced in Fig. |^ b) line-outs of the projection of the true density map n^^"*^ 
and of n"*" and n~ at x = lO/im from the simulated target. As it is evident, not only none 
of them reasonably reproduces the true density contour but their shapes also differ very 
much. As a result, standard Abel Inversion is not applicable in this case in order to produce 
a (reasonably) accurate density map. In Fig. ^ a) the line-outs of the projection of the 
generalized Abel inversion n^'^"''^'^ and the standard inversion of 50^ are confronted with the 
true density contour (dotted line) and uqi = uoPq + riiPi (dashed line). Standard Abel 
inversion applied to the mean phase-shift map still fails in reproducing a reasonable density 
map, while generalized Abel inversion gives us a contour which is everywhere well overlapped 
to the true one. Now, as a result, we can claim that with the use of the generalized Abel 
Inversion a good estimation of the simulated density map is achieved. 

Now, let us test the new algorithm in a physical condition which is often experimentally 
observed: a density with an axially symmetric background to which it is added an axially 
symmetric channel whose symmetry axis is not aligned with the one of the background (see 
Fig. ^ a)). As before, the maximum electronic density is well below the critical density 
^mctaf = 0.1 and the probe wavelength is l/zm. 

As in the previous example we apply the standard and the generalized Abel inversions to 
the phase-shift maps extracted by the interferogram reproduced in Fig. § b). In Figg. ^ c) 
and d), the best output of the standard Abel inversion (the one obtained with 50+) and of 
the generalized Abel inversion are shown, respectively. As in the previous example, standard 
Abel inversion produces a poorly accurate density map. A more quantitative analysis can 
be performed with the help of line-outs reported in Fig. ^, which shows that standard Abel 
inversion results should be rejected. On the contrary, generalized Abel inversion produces 
reasonably good results in almost all the density map but a thin band near the best symmetry 
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axis, where the dependence of n^^^"^ on cos(^) is much more complex then hnear. 

We face now with a noisy phase map in order to compare the noise content in the standard 
and in the generahzed Abel inversions. Here we will focus only on Gaussian and uncorrelated 
(white) noise, which is added to the phase-shift map of the interferogram reproduced in Fig. 
^ b) (see Fig. ^ (a)). To visuahze the noise which is propagated to the density maps n^"^^^ 
and n^°^'^^ (see Fig. |^ (b)), we subtract them to the density maps no and ni we have 
previously computed with the phase map of interferogram in Fig. ^ b) (in which no noise 
were introduced). The resulting error maps Stiq = uq — n^"*"*^ and 5ni = ni — n^°^^^ are 
finally confronted. Since the isotropic component of the density map ra^"*'^^ coincides with 
the standard Abel inversion of the symmetrized phase map, we can to compare the noise 
content in the standard and in the generalized inversions by simply comparing the noise in 
^Noise g^j^^ ^Noise^ ^'ig. ^c) line outs of Shq and 5ni are reproduced. 

Since the density maps are obtained integrating the uncorrelated noise with a kernel 
1/ \/C,'^ — which is rising in approaching the symmetry axis (r — > 0), we expect a noise 
sequence with a stronger component near r = 0. Futhermore, in computing rii (see Eq. |0.8| ) 
we face with the derivative of go that the 1/C term will contribute to enhance the 

noise in the r ^ region. Nevertheless, because of the regularization induced by the overall 
multiplication by r, a balancing of the two effects occur and the noise observed onto the 
^Noise j-Qg^p jg comparable with the one of the standard Abel inversion n^"^'^^ map, as it is 
clear in Fig. ^(c). 

To conclude, the generalized Abel inversion method we propose is very simple and effec- 
tive, it uses consistently the information carried by the whole phase-shift map and, as shown 
in the examples, it can be successfully applied to asymmetric cases for which the standard 
method based on only half-space phase shift map, fails. 
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Figures Caption 



Figure 1: a) Formation of an interferogram. The symmetry axis is x and the phase-shift is 
obtained integrating over the y direction, b) Because of an integration along the y axis, no 
departure of rig from a mirror-symmetric distribution can be detected. 

Figure 2: a) The z—x plane projection of the simulated electronic density. The radial profile 
is of the form n^™^ = uoPq + n-i Pi +^2^2, while the density is exponentially decreasing 
in the longitudinal direction, b) The simulated interferogram obtained with the density 
map n'^^^^ and Ap = 1 fim. c) and d) Projections onto the z — x plane of the density maps 
obtained with the standard and generalized Abel inversion of the phase-shift of the simulated 
interferogram b). 

Figure 3: a) Line outs of the z — x projection of the true density profile (dotted line), 
the sum of the Pq and Pi terms (dashed line) and of density profiles obtained via standard 
Abel inversion of (the mean of Scf)'^ and and of the generalized Abel inversion. The 
generalized Abel inversion is considerably more accurate than standard inversion, b) Line 
outs of the z — x projection of the density profiles obtained via standard Abel inversion of 
50"*" and S4>~. None of them well reproduces the true density profile (dotted line). 
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Figure 4: a) The z — x plane projection of the simulated electronic density. The radial 
profile is the sum of a background and a channel not aligned with its symmetry axis, while 
the density is exponentially decreasing in the longitudinal direction, b) The simulated in- 
terferogram obtained with the density map 'nF''"^^ and Ap = 1 ixm. c) Projection onto the 
z — X plane of the density maps obtained with the standard Abel inversion of 50"^ and d) 
projection of generalized Abel inversion. 

Figure 5: a) Line outs of the z — x projection of the true density profile (dotted line) 
and of density profiles obtained via standard Abel inversion of 5(j)s (dashed line) and of the 
generalized Abel inversion (continuous line). The generalized Abel inversion is considerably 
more accurate than standard inversion, b) Line outs of the z — x projection of the density 
profiles obtained via standard Abel inversion of 50+ (dashed fine) and 50" (continuous line). 
None of them well reproduces the true density profile (dotted fine) . 
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Figure 6: a) Line out of the noisy phase-shift map, obtained by summing up the phase-shift 
map of the interferogram in Fig. ^ and a Gaussian white noise map. b) Line outs of the 
resulting density maps n^"^'^^ and n^°^'^^ obtained via generahzed Abel inversion. Line outs 
of the hq and rii maps generated by inverting the noise free map are reported as a reference, 
c) Line outs of Suq and 5ni. The rms of the two noise sequences is comparable so the noise 
content of 71^°^^^ and ni°'^'"^ is similar. 
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